First Evidence of a Precessing Jet Excavating a Protostellar 

Envelope 

Jason E. Ybarra-*^ & Mary Barsony^'^ 

Department of Physics & Astronomy, San Francisco State University 
1600 Holloway Drive, San Francisco, CA 94132 

jybarraSstars . sfsu.edu, mbarsonySstars . sfsu.edu 

Karl E. Haisch, Jr.^ 

Utah Valley State College, 800 West University Parkway, Orem, UT 84058 

haischka@uvsc . edu 

Thomas H. Jarrett^ 

Infrared Processing and Analysis Center, MS 100-22 Caltech, Pasadena, CA 91125 

j arrett@ipac . caltech . edu 
Raghvendra Sahai^ 

Jet Propulsion Laboratory, MS 183-900, 4800 Oak Crove Drive, Pasadena, CA 91109 

Raghvendra . Sahai® j pi . nasa . gov 
and 

Alycia J. Weinberger 

Department of Terrestrial Magnetism, Carnegie Institution of Washington, Washington, 

DC 20015 

weinbergerOdtm . ciw . edu 
ABSTRACT 

^Observations with the Palomar 5 m telescope were obtained under a collaborative agreement between 
Palomar Observatory and the Jet Propulsion Laboratory. 

^and Space Science Insthute, 4750 Walnut Street, Suite 205, Boulder, CO 80301 



- 2 - 



We present new, sensitive, near-infrared images of the Class I protostar, Elias 
29, in the Ophiuchus cloud core. To explore the relationship between the infall 
envelope and the outflow, narrowband H2 1—0 S(l), Br7, and i^cont filters were 
used to image the source with the Wide-Field Infrared Camera on the Hale 5.0- 
meter telescope and with Persson's Auxiliary Nasmyth Infrared Camera on the 
Baade 6.5-meter telescope. The source appears as a bipolar, scattered light neb- 
ula, with a wide opening angle in all filters, as is typical for late-stage protostars. 
However, the pure H2 emission-line images point to the presence of a heretofore 
undetected precessing jet. It is argued that high velocity, narrow, precessing jets 
provide the mechanism for creating the observed wide- angled outflow cavity in 
this source. 

Subject headings: ISM: jets and outflows — stars:individual(Elias 29) — stars:formation 

1. Introduction 

Powerful, bipolar supersonic outflows are a hallmark of the earliest protostellar stages. 
Molecular outflows generally exhibit wide opening angles and velocities of ~ 10-30 km s~^. 
Optical jets powered by protostars exhibit small opening angles and velocities of ~ 100-400 
km s~^ (Mundt et al. 1990; Fridlund et al. 2005). Due to the disparate spatial and velocity 
distributions of the molecular and optical outflow components, a two-wind structure has 
often been invoked to explain the instances in which both outflow components are observed 
to be powered by a single source (Stocke et al. 1988; Davis et al. 2002; Liseau et al. 2005). 
However, an alternative scenario, in which optical jets and molecular outflows result from a 
single, coUimated flow is gaining wide acceptance (Hartigan et al. 1994; Sandell et al. 1999; 
Wolf-Chase et al. 2000; Yun et al. 2001). 

An outstanding, unsolved problem with jet-driven outflow models is how jets produce the 
wide-angle outflow cavities observed in both protostellar envelopes and in molecular outflows. 
A proposed solution is that of jets changing direction (Masson & Chernin 1993), for which 
there is now growing evidence (Bence et al. 1996; Reipurth et al. 1997; Barsony et al. 1998; 
Arce & Sargent 2005). Nevertheless, the spatial relationship between the outflow cavity and 
the currently active wandering jet has not yet been directly observed within the boundaries 
of a protostellar envelope. We use sensitive, narrowband near-infrared imaging to study 
this relationship: the scattered light emission from the protostellar envelope is imaged with 
a narrowband continuum filter, whereas the currently active jet regions within the infall 
envelope are traced in continuum-subtracted, narrowband H2 images, sensitive to shock- 
excited emission. Based on such observations of the Class I protostar, Elias 29=WL15 
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(q;2ooo = 16'*27"*09.43^ ^2000 = -24°37'18'/7) in the nearby (d=125 pc) p Ophiuchi molecular 
cloud, we report the first detection of a precessing jet carving out a protostellar envelope's 
cavity. 

2. Observations and Data Reduction 

Near-infrared, narrowband imaging of Elias 29 was undertaken at two facilities, Palomar 
and Las Campanas Observatories. The Wide Field Infrared Camera (WIRC) on the Hale 
5.0-meter telescope at Palomar Observatory, with its 8f7 x 8!7 field of view, was used to 
image the field including Elias 29 (Wilson et al. 2003). WIRC employs a 2048x2048 pixel 
Hawaii-II detector with a plate-scale of 0'.'2487 pixel"^. The K-band seeing was ~ I'.'O on 
the night of the observations, 2004 July 11. Two narrowband filters were used for imaging: 
an H2 1-0 S(l) filter centered at 2.120 //m with a 1.5% bandpass and a line-free continuum 
filter (Kcont) centered at 2.270 //m with a 1.7% bandpass. The telescope was dithered in a 
nine-point pattern, with 6" offsets between dither positions. The integration time at each 
dither position was 60 sec (30 sec x 2 coadds). The dither pattern was repeated, after a 5" 
offset of the telescope, until a total integration time of 27 minutes was reached in each filter. 
Elias 29 was rc-observed with Persson's Auxiliary Nasmyth Infrared Camera (PANIC) on 
the Baadc 6.5-meter telescope at Las Campanas Observatory (LCO) on the night of 2005 
June 17 UT. The K-band seeing was 0'.'5. PANIC employs a 1024x1024 Hawaii-I detector 
with a plate scale of 0'.'125 pixel^^ (Martini et al. 2004). The two narrowband filters used for 
imaging at LCO were the H2 1-0 S(l) filter centered at 2.125 //m with a 1.1% bandpass, and 
the Br7 filter, centered at a wavelength of 2.165 jim with a 1.0% bandpass. The Br7 filter 
served as a surrogate narrowband continuum filter for the PANIC dataset, in the absence 
of a fine- free, narrowband continuum filter adjacent to the H2 filter at LCO. The observing 
pattern for Elias 29 at LCO consisted of interleaved on- and off-source integrations, of 30 sec 
duration each, through each filter. On-source integrations were dithered in a checkerboard 
pattern, with 12" offsets. Five-point dithered sky observations were taken in a clear patch 
3' E of Elias 29. The final on-source integration times with PANIC were 70 mins. through 
the H2 filter and 67.5 mins. through the Br7 filter. 

All data were reduced using IRAF'^. Object frames were sky-subtracted, flat-ficldcd, 
and corrected for bad pixels for each night's observations through each narrowband filter. 
Individual processed images for each filter/instrument combination were then aligned using 
the available point sources in common with the IRAF task IMCENTROID, before combining 
to produce final images. 

Although the H2 filters are so-called because they transmit the 2.12/im emission line. 
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they both also transmit some continuum emission, if present in the source. Therefore, to 
to trace the morphology of any outflow/jet component within the protostellar envelope of 
Elias 29, the bright continuum emission from the scattering envelope must be subtracted 
from the H2 filter images. This requires proper scaling and alignment of the narrowband 
"continuum" images before subtracting them from the respective H2 filter images. These 
tasks were achieved separately for the WIRC and PANIC data. For the PANIC images, 
shown in Figure 1, the Br7 filter image was scaled by a factor of 0.82, which is simply the 
ratio of the transmission of the H2 filter to that of the Br7 filter, from the filter transmission 
curves. Proper alignment of the H2 and scaled "continuum" images before subtraction, 
crucial for obtaining accurate pure H2 line-emission images, was achieved by using the IRAF 
task, IMCENTROID, for common point sources in each filter. 

For astrometry, a plate solution was determined for the WIRC data by using a least 
squares fit for 7 point-like objects in common with 2MASS ^. The 2MASS reference objects 
had uncertainty ellipses of ~ Of! 09. The plate solution had a Of! 2 RMS residual, our estimated 
astrometric uncertainty. 

3. Results and Discussion 

The reduced WIRC and PANIC images of Elias 29 appear essentially identical. Only 
the PANIC images are presented here, since these were obtained in better seeing and have 
higher signal-to-noise. Figure 1 shows our results, all presented at a linear greyscale stretch 
to emphasize the presence of the pure emission features, 3a and 3b in Figure la, and their 
absence from Figure lb. Figure la shows the appearance of Elias 29 through the narrowband 
H2 filter. Figure lb shows Elias 29 through the surrogate "narrowband continuum" (Br 
7) filter. Since the NIR continuum emission overwhelms the strength of the pure H2 line 
emission for this object, Elias 29 looks similar in both filters. Figure Ic shows the continuum- 
subtracted, pure H2 emission-fine image of Elias 29. Saturation effects from the extreme 
brightness of Elias 29 prevent correct image subtraction within the central ~ 5''5 radius 
region, resulting in the black and white circular artifacts at the center of Figure Ic. Outside 
this area, five regions of pure H2 line emission are discovered, having been missed by previous 
investigators (Gomez et al. 2003; Khanzadyan et al. 2004). The appearance of Elias 29 in 



■^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

^http://www.phys.vt.cdu/~jhs/SIPbeta/astrometrycalc.html 
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Figures la & lb is consistent with scattered light models of Class I protostars consisting of an 
infall envelope containing a bipolar cavity (Whitney et al. 2003). The boundary of the wide- 
angled cavity walls within the protostellar infall envelope is indicated by the lowest contour 
levels overplotted on Figures lb & Ic. The lack of a narrow "waist" or distinctive hourglass 
shape at this angular resolution is consistent with a disk of small spatial extent (B. Whitney, 
priv. comm.). Elias 29 appears strikingly different in the continuum-subtracted, pure H2 
emission line image than in the narrowband filter images. The pure H2 emission objects, 
labelled 1, 2a, 2b, 3a, and 3b, appear at the same location, and have similar morphologies 
in both the PANIC and the WIRC (not shown here) continuum-subtracted images. Taken 
together, the H2 emission knots form a sinuous structure with S-shaped point-symmetry 
about the center. 

Table 1 hsts the observed characteristics of each H2 feature. The first column lists the 
feature designations from Figure Ic. The second and third columns list the coordinates of 
each H2 emission region, corresponding to the WIRC pixel with the highest counts in each 
object. The coordinates are derived from the WIRC data, since there were not enough 
sources in the smaller PANIC field to provide a plate solution for absolute astrometry. The 
final, continuum-subtracted PANIC image was used to determine the signal to noise ratio 
of each emission feature, listed in the last column of Table 1, since the seeing and signal 
to noise were better in the PANIC than in the WIRC dataset. A rectangular aperture (of 
dimensions given in the fourth column of Table 1) was used to determine the average counts 
in each feature. Several background regions were used to calculate the average background 
noise level and its standard deviation. The backgound level was subtracted from the average 
counts in the rectangular aperture containing each II2 emission region, and the result was 
divided by the standard deviation of the background to give the SNR. Only features 3a and 
3b have cr < 3. Nevertheless, these objects arc real, since they appear at identical locations 
and have similar morphologies in both datasets (WIRC and PANIC). 

H2 emission associated with molecular outflows generally arises from shock interactions 
of the unseen, fast jet/wind component with the ambient medium. Elias 29 is known to 
drive a CO outflow (Bontemps ct al. 1996; Sekimoto et al. 1997), with velocities as high 
as ~ 80 km (Boogert et al. 2002), although its morphology is unknown at the scales 
of the images presented here. In some Class I protostellar envelopes, both spatially dif- 
fuse H2 emission and compact H2 emission knots have been identified in the same source 
(Yamashita & Tamura 1992; Davis et al. 2002). In such cases, the diffuse component is 
found to trace the infrared continuum reflection nebula, whereas knotty, compact emission 
traces the fast jets observed at optical or cm-continuum wavelengths. By analogy, the knotty 
H2 emission features identifled in Figure Ic are likely to be associated with the jets driven 
by Elias 29. Since the cooling time of shocked 2.12 /im H2 emission is of order ~ a few years 
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(ShuU & Beckwith 1982; Allen & Burton 1993), the newly identified features in Figure Ic 
arc tracing recent shock activity. The three H2 emission features in Figure Ic: 1, 2a, and 
2b, describe an S-shaped, point-symmetry, about the center, consistent with a precessing 
jet. Features 3a and 3b are extended and diffuse, close to the inferred cavity wall within 
the southern hemisphere of the protostellar envelope of Elias 29. More quantitative results 
on the precessing jet, such as its precession period and initial opening angle await HST 
NICMOS polarimetry. Ground-based, AO-assisted integral field spectroscopy of the H2 hne 
would establish the velocity field. 

Hydrodynamic simulations of molecular outflows driven by pulsed, precessing protostel- 
lar jets have become available only recently (Lim 2001; Rosen & Smith 2004; Smith & Rosen 2005). 
One stated aim of these studies is to model the outflows to infer properties of the (unseen) 
driving jets, such as jet power and precession. Currently, only a limited set of models have 
been calculated with accompanying simulated H2 images. Assumed jet radii are 1.7 x 10^^ 
cm, corresponding to ~ I'.'O at the distance to Elias 29 (so, well beyond the disk region where 
jets may be generated). The pulsation is input as a 30% velocity modulation with a 60 yr 
period into the models. The timespan over which the models are calculated is short (t < 
500 yr), but much longer than the dynamical time of a few decades for a jet travelling a 
conservative ~ 100 km s~^ to reach the H2 Features 1 and 2a in Elias 29. This, compared 
with the very fast coohng time of H2, validates comparison of the Ehas 29 pure emission-line 
H2 image (Figure Ic) with simulated H2 images. 

Models of molecular jets with precession angles of 5°, 10°, and 20° have been calcu- 
lated, although the precession angle of the Elias 29 jet may well be even wider, based on 
the projected half-angles formed by Features 1 and 2a with the central symmetry axis of 
the observed reflection nebulosity. Simulated H2 images from two general classes of pre- 
cessing jet models have been published: Jets with assumed precession periods of 50 yr, 
termed "fast" precessing, and jets with precession periods of 400 yr, termed "slow" pre- 
cessing (Rosen & Smith 2004; Smith & Rosen 2005). The terms "fast" and "slow" here are 
relative to the flow evolution time of 500 yr. Comparison of Figure Ic with synthetic H2 
f = 1 — S(l) model images from such simulations leads to the conclusion that this jet 
exhibits slow, rather than fast, precession, that is to say, qualitatively. Figure Ic resembles 
the H2 1^0 morphologies presented in Figures 6 & 7, rather than those of Figure 5, of Smith 
& Rosen (2005). For slow precession, the simulations show ordered chains of bow shocks and 
meandering streamers of H2, which contrast with the chaotic H2 structures produced by jets 
in rapid precession. The production of specific simulated H2 images for direct comparison 
with Figure Ic is beyond the scope of this paper, but may be attempted in the future. 

Images of both scattered hght cavities and H2 jet emission within them exist for only 
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a handful of Class I objects (Davis et al. 2002). For all of these, the jets he well within 

the outflow cavity walls, close to, if not precisely along, the symmetry axis of each outflow 
cavity. By contrast, the jet of Ehas 29, as traced by H2, is found close to, and, in one case 
(Feature 2b) along, the cavity walls (Fig. Ic). 

Examples of other S-shaped, precessing outflows are known (Gueth et al. 1996; Schwartz 

& Greene 1999; Hodapp et al. 2005; Sahai et al. 2005). However, no precessing pro- 
tostcllar jets have been imaged at their base within their scattered light cavities, either 
because they are too embedded to be detected in the NIR (LI 157, HH211-mm), or be- 
cause of lack of narrowband continuum data excluding the H2 emission line for proper im- 
age subtraction(IRAS 03256+3055). Binarity has inevitably been invoked to account for the 
presence of precessing jets (Terquem et al. 1999). However, Elias 29 is known to be single 
(Simon et al. 1995; Ghez et al. 1993). At scales between ~ 1 — 60AU, lunar occultation 
observations show that 90% of the K-band flux from Elias 29 originates in a component ~ 7 
mas diameter (~ 0.8 AU) with the remaining 10% of the K-band flux originating in a much 
larger, 415 mas diameter, diffuse component (Simon et al. 1987). It is then inferred that 
within 0.8 AU, Elias 29 is a single object. Alternative scenarios leading to jet precession 
that may be operative in Elias 29 might be disk oscillations associated with a propeller-driven 
outflow, recently discovered in MHD simulations (Romanova et al. 2006) or warping of the 
accretion disk caused by the magnetically driven jets (Lai 2003). Future high- resolution 
integral field spectroscopic observations of both the jet and the disk components could test 
the latter model, which predicts the disk and jets to precess in opposite senses. 
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NSF AST-0206146 to M.B. and by the NASA LTSA Program (NAG5-8933). NASA/LTSA 
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Table 1. H2 emission features 



H2 Feature 


Q!2000 


(^2000 


Aperture 


S/N 


Designation 


(h m s ) 


/ // 






1 


16 27 08.79 


-24 37 16.5 


0'.'6 X 0'.'5 


3.5 


2a 


16 27 10.04 


-24 37 20.8 


0'.'6 X 0'/5 


3.9 


2b 


16 27 10.70 


-24 37 24.3 


1'.'3 X 0'/5 


14.2 


3a 


16 27 11.64 


-24 37 56.1 


0'.'6 X 0'.'5 


2.6 


3b 


16 27 11.88 


-24 38 00.0 


O'.'g X 0'.'4 


2.8 
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Fig. 1. — Narrowband Images of Elias 29. N is up and E is left in all panels. Offsets from 
RA(2000)= 16'*27"'09.43^ and Dec(2000)= -24°3ri8'.'7 are indicated, in arcsec units, (a) 
Elias 29, observed through the narrowband H2 filter with PANIC, (b) Elias 29 imaged with 
the Br7 filter with PANIC, contour levels are 3, 5, 10, 20, 50, 500 a. Note the absence of 
Features 3a and 3b in the Br7 image, and their presence in the II2 filter image, (c) Pure 
H2 emission-line image of Elias 29, obtained by subtracting the scaled and aligned Br7 filter 
image from the PANIC H2 filter image with the Br7 10 cr contour superimposed. The pure 
H2 emission line features discussed in the text and hsted in Table 1 are labelled. 



